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Theory of Liquid Crystal Anchoring at a Porous Surface 

David L. Cheungi'U 

'Department of Physics and Centre for Scientific Computing, 
University of Warwick, Coventry, CV4 7AI, UK 

Using classical density functional theory (DFT) the effect of bringing a liquid crystal (LC) into 
contact with a porous substrate or matrix is investigated. The DFT used is a combination of the 
Onsager approximation to evaluate the excess free energy of the LC fluid and quenched annealed 
DFT to evaluate the interaction between the fluid and the porous substrate. When the fluid align- 
ment far from the substrate is held perpendicular to its surface there is a thin layer of fluid aligned 
parallel to the substrate surface for low matrix densities. For higher matrix densities this director 
deformation propagates into the bulk of the fluid. Consideration of a system without conflning 
walls suggests that for low matrix densities normal alignment is metastable compared to parallel 
alignment, while for higher matrix densities it is unstable. 



I. INTRODUCTION 

The interaction between liquid crystals (LC) and 
solid substrates of of great interest, both scientific and 
technologicali. The presence of the surface both breaks 
the symmetry of the LC phases and often leads to align- 
ment in a given direction. This tendency, often called 
anchoring, is vital to the use of LC in display applica- 
tions and can be profoundly influenced by the structure 
of the surface. 

Commonly surfaces may be rough or porous, e.g. 
Si02i^, which can lead to deviations from the behaviour 
expected from smooth surfaces^. A porous substrate or 
matrix may be regarded as a system of immobilised par- 
ticles. Recent simulations of LC near rough walls^ using 
such a model have shown that properties such as the an- 



choring coefflcient and transition pressures of a confined 
LC are influenced by the roughness of the substrate. Due 
to the computational expense only a few surfaces at a 
given roughness were studied. It would thus be desirable 
to study a larger (ideally infinite) number of surfaces. 

One route to this is through replica or quenched an- 
nealed density functional theory (QA-DFT)i2i24i2i. In 
this the system comprises two components: the first, the 
quenched component, models the porous substrate, the 
second, the annealed component, models the fiuid. The 
density distribution of the matrix, averaged over all ma- 
trix realisations, is determined though minimisation of a 
grand potential in the absence of the fiuid. The density 
distribution of the fiuid is then found by minimising a 
grand potential containing both matrix and fiuid densi- 
ties. This theory is fully outlined in the following section, 



along with details of the systems studied. The results 
that ensue from applying it to LC anchoring at porous 
surfaces are outlined in Sec. IIIII Finally brief summary 
of this work and suggestions for future work are given in 
Sec.HVI 

II. THEORY AND MODEL 
A. Density functional theory 



the fluid and matrix together, i.e. ^[pynif, u), pf{r, u)] = 
^miPmir , u)] + Q{[pai{r , u) , p{{r , u)]"^ . Here we are in- 
terested in uniaxial molecules characterised by a po- 
sition r and orientation u. For an ideal matrix (i.e. 
with vanishing interactions between the matrix particles) 
ilm[pm{r,u)] is 



In QA-DFT the grand potential is the sum of the grand 
potential of the matrix alone and the grand potential of 
I 



f3n, 



,[pmir,u)] = / drdup^{r,u) {log pir,{r,u) - l + /3V;„(r,M) - (3p.m} 



(1) 



r 



where V„i{r,u) is the external potential acting on the 
matrix, /im is its chemical potential, and f3 = l/k^T. If 
the external potential Kn(r,M) is used to confined the 
matrix to a region of space then the matrix grand po- 



tential may be minimised analytically to give Pni{r, u) 
Pni = exp(/3/iin) inside this region and outside it. 

The grand potential for the fluid component is 



f3n{[pf{r,u),p^{r,u)] = / drdupf{r,u) [log pi{r, u) - I] + P drdu [V {r , u) - p] pf{r , u) 



+ I3F,^ [Pf(r, u)] -f PF^J [pf(r, u), p^{r, u)] 



r 



(2) 



where pf(T", u) is the density distribution of the fluid com- 
ponent, V{r,u) is the external potential acting upon it, 
and p is the fluid chemical potential. F^^ [pf(r,M)] and 
F^^ [pi{r, u), pmir, u)] are the excess free energies due to 



interactions between the fluid molecules and fluid-matrix 
molecules respectively. From hereon in the subscript f de- 
noting the fluid component will be omitted. The excess 



FE is evaluated within the Onsager approximationii 

(3Fc^[p{r,u)\ = -- driduidr2du2 

f{ri2,Ui,U2)p{ri,Ui)p{r2,U2) (3) 

where /(ri2, 111,1*2) is the Mayer function. For conve- 
nience we assume that the interaction between the matrix 
and fluid particles is identical to the interactions between 
the fluid particles. 

In this work we are only concerned with systems that 
vary in the z direction only; the surfaces are assumed to 
be homogeneous in the x and y directions. The method 
I 



for finding the equilibrium density is the same as in previ- 
ous workiSiii; the density and its logarithm are expanded 
in spherical harmonics 



p{z, u)^Y^ pin-,{z)Y;„^{u) (4a) 

Im 

log p{z,u)^^Pi,n{z)Yini{u) (4b) 

Im 

where Yfm(M) is a spherical harmonic. The Mayer func- 
tion is expanded in rotational invariants. Inserting these 
into the grand potential (Eq. ^ and integrating over an- 
gles and the x and y axes gives 



f3Qf[p{r,u),p^{r,u)] = / dz ^pimiz) lpem{z) + pVimiz) - V^{1 + (3p)Seo\ 

- / dzidz2 y^ Ai^i^rn{zi~ Z2)pi^m{zi)\pi2m{z2) j=Pm{z2)5l>^Q 



ilt^ra 



(5) 



where Vim{z) are the spherical harmonics coefficients of 
the external potential and Ai-^i^^lzi — Z2) are the spheri- 
cal harmonics coefficients of the excluded area^^. Pni{z) is 
the matrix density, equal to p^ for Zmin < z < z,nax and 
otherwise. The equilibrium fluid density is then found by 
numerically minimising this with respect to pi-ai{z) using 
the conjugate gradients methocU^. When needed pira{z) 
are found from Eqs. ^ Once the pira{z) that minimises 
fif has been found the fluid density profile is given by 
p{z) — J du p{z,u). The orientational ordering of the 



fluid is described through the ordering tensor 



Qal3{z) 



du p{z,u)uaUi3-'^Sai3, a,P = x,y,z . 

(6) 



The order parameter is given by the largest eigenvalue 
of Qapiz) and the director by the eigenvector associated 
with it. 

Both the fluid and the matrix particles are modelled as 
hard ellipsoids of revolution of elongation e = a/h = 15. 
The chemical potential is set to // = 2.0kB well inside 
the nematic phase. This class of model have been well 
studied as model liquid crystals^^. Two different systems 



were considered. In the first the fluid is confined between 
two hard walls at z = and z = L = 2005, with the 
matrix confined in the region < z < Zmax — 306. The 
wall at z = L is used to provided (homeotropic or planar) 
alignment far from the matrix. In the second system the 
matrix is placed in the centre of a periodic system of 
width 2005. The fluid density is then minimised from 
starting states parallel and normal to the matrix surface. 

III. RESULTS 
A. Confined Geometry 

Shown in Fig. ^ are the density and order parame- 
ter profiles for the slab geometry with both planar and 
homeotropic alignment at the far wall. In all cases far 
from the matrix p{z) is constant. For the lowest ma- 
trix density (p,n = 0.055^'^) the fluid density with the 
matrix is non-zero, although far lower than in the fluid 
outside the matrix. At larger pm the fluid density inside 
the matrix is zero. The variation in p{z) near the matrix 
surface depends on p,„ and the alignment far from the 
matrix. For planar alignment (Fig.^) and p,n < 0.255"'^ 
p{z) drops smoothly from the bulk value at z « 455 (one 
molecular length from the matrix surface) to just inside 
the matrix. For larger values of pm a peak appears in p{z) 
at z « 37.56 (half a molecular length from the surface), 
corresponding to an absorbed layer of molecules. p{z) 
then drops sharply to about one molecular diameter 
outside the matrix. For these pm the matrix becomes im- 



penetrable to the fluid and starts to behave like a rough 
hard wall*. When the far wall gives rise to homeotropic 
alignment p[z) has a slight peak at z « 406 for all pm 
that grows stronger for increasing p^. For the highest 
values of pi„ studied however, this peak is weaker than 
in the corresponding system with planar alignment. 

Shown in Fig. ^ are order parameter profiles (>S'(z)) 
for the confined LC with a far planar wall. In common 
with the density profiles these are constant far from the 
matrix, while for all pm S{z) = inside the matrix. 
For Pm = 0.055^'^ S{z) drops smoothly to while for 
a slightly higher pm a sharp peak appears just inside the 
matrix surface. At higher pm S{z) also drops smoothly 
from its bulk value to 0, although this occurs outside the 
matrix and over a smaller range of z. S{z) for the far 
homeotropic wall is shown in Fig. ^. The behaviour of 
S{z) near the matrix surface differs markedly for different 
Pni. For low Pm there are two peaks in S{z) either side of 
the matrix surface, while deeper in the matrix S{z) goes 
to zero. At larger pm the order parameter drops to at 
the matrix surface. 

For the system with far homeotropic alignment the be- 
haviour of S{z) near the matrix surface in infiucnced by 
the behaviour of the director. Plotted in Fig. El is 0{z), 
the angle between the director and the z axis (the wall 
normal). For low pm the director lies parallel to the z- 
axis {0{z) = 0) far from the matrix surface. Near the 
matrix surface however, it abruptly changes to lie in the 
xy plane {0{z) — 90°). This preference for parallel align- 



ment may be understood as the matrix consists of disor- 
dered ellipsoids, so the interface between it and the fluid 
is equivalent to an isotropic-neniatic interface. Previous 
studies of this and similar modelsiStii have shown that 
the director prefers to lie in the plane of the interface 
For Pin > 0.26 the preference for planar anchoring at 
the matrix surface leads to continuous director variation 
through the cell from homcotropic at the far wall to pla- 
nar at the matrix surface. By contrast for the far planar 
wall is constant 9{z) = 90° for all pm- 



,18 



/37 



pn[pir,u)]-pnb[pbiu) 

2A 



(7) 



where fti, [/Ofc(l)] is the bulk grand potential. This is plot- 
ted in Fig. ^ for both initial parallel and normal align- 
ments. As may be seen for pm < 0.256^'^ 7 is lower for 
parallel alignment, indicating that normal alignment is 
metastable. For higher pm normal alignment is unstable 
and 7 is identical for both initial conditions. 



B. Open Geometry 

Density and order parameters for the open system are 
shown in Fig. El As can be seen for both the homeotropic 
and planar anchoring p(z) and S{z) are similar to those 
for the confined system: far from the matrix these are 
constant while close to the surface of the matrix these 
drop to smaller values. Also, with the exception of the 
Pm = 0.056^^ system, p(z) drops to within the matrix. 

The behaviour of the director is different in the open 
system. When the initial configuration for the minimisa- 
tion has the director aligned parallel to the matrix surface 
the director remains in this state. For an initial normal 
alignment the final state also has normal alignment when 
Pm < 0.25&~^, while for higher pi„ the final state has par- 
allel alignment. This suggests that the normal alignment 
is metastable with respect to parallel alignment. In or- 
der to confirm this we calculate the surface free energy 7, 
given by the excess (over bulk) grand potential per unit 



IV. CONCLUSION 

In this paper the effect of bringing a liquid crystalline 
fiuid into contact with a porous substrate has been stud- 
ied using QA-DFT. It was found that aligning the LC 
parallel to the substrate surface is favourable to normal 
alignment. For matrices of low density normal alignment 
is found to be metastable, while for higher pm it is un- 
stable. This preference for parallel alignment is in agree- 
ment with studies of the nematic-isotropic interface for 
the same model. 

There are a number of possible extensions to this work. 
It would be of interest to compare the QA-DFT results 
to those of computer simulations. Also interesting phe- 
nomenon in the vicinity of the nematic-isotropic transi- 
tion that may be examined or for LC mixtures and for 
different interactions between the fluid and matrix par- 
ticles. 
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Figure Captions 

Fig. 1. (a) Density profiles {p{z)) for confined LC with 
far planar wall. Solid line pi„ — 0.056^^, dot- 
ted hne p,„ — 0.106^^, dashed hne pm = 0.256^^, 
long dashed line pm = 0.30&^'^, and dot-dashed 
line pni — 0.406^^. (b) Density profiles for con- 
fined LC with far homeotropic wall. Symbols as in 
(a), (c) Order parameter profiles {S{z)) for con- 
fined LC with far planar wall. Symbols as in (a), 
(d) Order parameter profiles for confined LC with 
far homeotropic wall. Symbols as in (a). 

Fig. 2. Director angle 9{z) for confined LC with far 
homeotropic wall. Solid line shows pm = 0.056^'^, 
dotted line shows pm = 0.256^"^, dashed line shows 
Pni — 0.266^^, long dashed line shows pi„ — 0.306^^ 
and dot-dashed line shows pm = 0.40&~'^. 

Fig. 3. (a) Density profiles for open system with parallel 
alignment. Solid line pni = 0.056^^, dotted line 
Pni = 0.10&^^, dashed line p^ = 0.256"^, long 
dashed line pm — 0.306^'^, and dot-dashed line 



Pm = 0.406~^. (b) Density profiles for open system 
with perpendicular alignment, (c) Order parame- 
ter profiles for open system with parallel alignment. 
Symbols as in (a) . (d) Order parameter profiles for 
open system with perpendicular alignment. 

Fig. 4. Surface free energy (7 = Q — Q,b) per unit length for 
open systems with parallel (solid line) and normal 
(dashed hue). 
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